Net Charge Fluctuations as a signal of QGP from 
Polyakov— Nambu— Jona-Lasinio model 
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We report the first model study of the net charge fluctuations in terms of D — measure within 
the framework of the Polyakov— Nambu— Jona-Lasinio model. Net charge fluctuation is estimated 
from the charge susceptibility evaluated using PNJL model. A parameterization of the freeze-out 
curve has been used to obtain D as a function of <Js. We have discussed our results vis-a-vis recent 
experimental findings from ALICE collaboration. 
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Strongly interacting matter at very high temperatures 
and densities is expected to undergo a transition from 
confined state of colored charges, the hadronic phase with 
broken chiral symmetry to a partonic phase in which chi- 
ral symmetry is restored and/or quarks are deconfincd 
A good understanding of this transition is relevant 
for studies in the fundamental interactions in particle 
physics, as well as for the physics of early universe and 
neutron stars. Thus it has become an issue of great in- 
terest in recent years, both theoretically and experimen- 
tally To this end, it is essential to identify unam- 
biguous signals which would establish the formation of a 
quark-gluon plasma (QGP). One such viable signal is the 
fluctuations of net electric charge Q 0, H| . It has been 
argued that this fluctuation is proportional to the square 
of the electric charge which takes up distinct values for 
the hadronic and QGP phases. While the unit of Q in 
the hadronic phase is 1, that in the QGP phase is 1/3. 
This may result in the fluctuation in net charge to vary 
with the change of phase, with the net charge remaining 
unaffected. Such fluctuations in heavy ion collision ex- 
periments are measurable via event-by-event (EbE) anal- 
ysis where one single event corresponds to a set of 
innumerable particles produced in a single collision of rel- 
ativistic nuclei. This method deals with measurement of 
a given observable on an EbE basis and study of fluctu- 
ations over an ensemble of events. 

Measuring charge fluctuations: To reduce systematic 
uncertainties of measurable quantities in heavy-ion ex- 
periments, it is useful to consider suitable ratios of quan- 
tities that are expected to have similar systematic be- 
havior. Here for measuring charge fluctuations a suitable 
observable could be the ratio 
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of net charge Q to total charge N c h- One could also 
consider a ratio R defined as, 
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Here, if one uses the approximation (Nch) » (Q), then 
(SR 2 ) ~ A(5F 2 ) where, 
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So the signal in fluctuations of F is amplified four times 
in the fluctuations of R. Now R is related to ratio of net 
charge fluctuation to (N c h) via a quantity D, defined as: 



D = (N ch )(SR 2 ) = 4 
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where n c h = (N c h) /V is the total charge density. n c h 
may be obtained by adding the contribution from parti- 
cle and anti-particle distributions while the dimensionlcss 
charge susceptibility \Q may be obtained from the pres- 
sure P of the system as, 
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where /xq is the electric charge chemical potential. 

This definition is useful over that of the ratio fluc- 
tuations as the general form of the latter may be 
quite complicated. As mentioned above that only if 
(N c h) >> (Q) the simplified relation exists as given 
above. In general this condition is expected to be sat- 
isfied only for a large T and very small \xb- Once the 
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fi b becomes large the net charge Q increases and the ap- 
proximation would fail. However the definition of D is 
quite general and holds even if not so simply related to 
the ratio fluctuation. 

A simple estimate of D was made considering the 
hadronic phase to be composed of pion gas and the par- 
tonic phase as computed in the Lattice Gauge Theory 
Q • This gave the value of D to be ~ 4 for the hadronic 
phase and ~ 1 for QGP phase. 

The first measurement of D in experiments have been 
reported recently by the ALICE collaboration ■ They 
have obtained the net charge fluctuations in a rapidity 
window 0.2 < An < 1.6 using center of mass energy y/s = 
2.76 TcV in Pb-Pb run for different centralities. The D 
for increasing An continues to fall and is just short of the 
saturation region as expected from UrQMD simulation 
results Q . The analysis of D and its variation with center 
of mass energies in the range of 19.6GcV < s/s < 200GcV 
has also been presented in the same report. The data 
used were obtained by the STAR collaboration . 

Results in PNJL model: Here we report on the study 
of net charge fluctuations in terms of the D — measure 
using the Polyakov loop enhanced Nambu— Jona-Lasinio 
(PNJL) model. In this model the quarks and Polyakov 
loop fields are the basic degrees of freedom. The quarks 
while interacting with the Polyakov loop also has a four- 
fermi self-interaction. Similarly the Polyakov loop fields 
interact via a Landau-type potential. The details of the 
model used for a 2 flavor system may be found in [l^, EH ■ 
The extension to 2+1 flavor system has been done in 
14j . Detailed studies of fluctuations and correlations of 
various conserved charges were performed with the PNJL 
model both for 2 flavor [H,[3 and for 2+1 flavor [TMl 
systems. 

To compute D we evaluate XQ and n ch using the PNJL 
model. The method of obtaining xq is quite standard 



as has been discussed by us earlier |15|. On the other 
hand n c h are to be calculated from the quark distribution 
functions as they appear in the PNJL model. 

The behavior of XQ and n c h/T 3 with T/T c are shown 
in FigQ] for various values of the baryon chemical po- 
tential us and for the cases of 2 flavor and 2+1 flavor 
systems. Here T c is the crossover temperature at the 
corresponding values of /ib- The quantities under con- 
sideration show qualitatively similar behavior. There is 
a sharp rise close to T/T c ~ 1 from a negligibly small 
value for low temperatures. This is followed by a satu- 
ration at high temperatures close to the corresponding 
values for massless free quarks. The values of XQ and 
n c h/T 3 , at any T/T c are larger for larger /xg. However 
as T increases they all tend towards the limit of free gas 
at ^b = 0. 

We now consider the quantity D/Df ree and study its 
temperature and density variations. Here Df ree (T, /is) is 
the temperature and chemical potential dependent limit 
of D(T,/ib) for a free massless gas of quarks. In Fig 121 
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FIG. 1: (color online) Variation of \Q an d n c h/T 3 (lower and 
upper set of curves) with T/T c for different values of fis for 
2 flavor (upper panel) and 2+1 flavor (lower panel). 



we show the variation of D/Df ree as a function of T/T c 
for both 2 and 2 + 1-flavor cases. We have chosen four 
representative values of /.is- It is observed that D al- 
ways remains above its free field limit and approaches 
this limit at high enough T. The sharpest transition oc- 
curs near T = T c . Now both xq an d n c h are smaller than 
their corresponding free field limit below T = T c for any 
/i b ■ Therefore D > D f ree implies that n c h is much more 
suppressed than xq m the confined phase. This effect 
seems to be much more prominent as [lb is increased. In 
passing, let us mention that (Q) / (N c h) in our studies 
varied from a value of 0.025 for high T and low /is to 0.3 
at the other extreme. 

Connection with heavy-ion collision experiments: It 
would thus be interesting to check what happens with 
further increase of [Ib- The variation of D/Df ree with 
fis at different temperatures arc shown in Fig. [3l Again 
we find D to remain above its free field limit for all T 
and hb- For a lower temperature ~ 100 MeV, there is 
an initial rise and a subsequent fall with increasing [is- 
This non-monotonic behavior at low temperatures may 
pose a problem in using D as a direct indicator of the 
phase of strongly interacting matter. 

Nonetheless, with an input of temperature and chem- 
ical potential from particle multiplicities at the freeze- 
out surface in heavy ion collision experiments, one may 
study the expected nature of D for different experimental 
conditions. The independent thermodynamic variables 
in the PNJL model are T, [is, f-Q, and jig, where the 
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FIG. 3: (color online) Variation of D/Df ree with \ib for three 
values of T around fis = for 2 flavors (upper panel) and 2+f 
flavors (lower panel) 



FIG. 2: (color online) Variation of D with T/T c for different 
values of [ib around (iq — 



latter is the strangeness chemical potential. Parameter- 
ization of the freeze-out conditions as a function of y/s 



are available in the literature (see e.g. [2(3, 21 1). For a 
given set of the thermodynamic variables we found the 
variations in y/s are within 10% for different parameter- 
izations. Here we choose the parameterization in [2l[ to 
model the freeze-out curve as: 



where B, Q and S indicates the values of d and e for the 
corresponding chemical potential given in Eqn(7l 




2 flavor ■ 
2+1 flavor ■ 



T(fi B ) = a-bn 2 B - cfi% (6) 

IJ>b,q,s(Vs) = — /= (7) 

1 + esfs 

where, a = 0.166 ±0.002 GcV, b = 0.139 ±0.016 GeV" 1 , 
c = 0.053 ± 0.021 GeV -3 and d and e are given as: 
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FIG. 4: (color online) D as a function of y/s computed along 
the freeze-out curve. 

We have thus calculated D (IS c\ function of y/s along 
the freeze-out curve. Results are shown in Fig. |U We 
have varied y/s from 5 GeV to 3 TeV. D picks up a value 
of ~ 4 for 2 flavor and ~ 3.5 for 2+1 flavor at the low 
y/s, drops down to a value of 2.6 for 2 flavor and 2.4 for 
2+1 flavor around y/s ~ 200 GeV, and saturates at these 
values even for increasing y/s. It is highly exciting to find 
that the general features of D vs y/s curve obtained in 
the PNJL model are found to be similar to those obtained 
directly in heavy-ion collision experiments as shown in 
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Fig. 4 of Ref . (lot . Furthermore the numerical range of D 
itself is exceptionally consistent. 

It should however be remembered that D as given in 
Fig. [~4]is the value obtained when the system is in com- 
plete thermal equilibrium at the given values of temper- 
ature and chemical potentials. Since here we are on the 
freeze-out curve, we are always inside the hadronic phase, 
i.e. the whole of the curve in Fig. |4] is corresponding to 
varying environmental conditions in the hadronic phase. 
Thus if our results become completely consistent with 
experimental results the outcome will be that there is no 
signature of partonic phase in D. At present it seems that 
the results for STAR data are above and those of ALICE 
data are below our model curve. A more concrete anal- 
ysis would be possible once the complete experimental 
data for D are published. 

To summarize, we presented here model study of net 
charge fluctuations in terms of .D-measure using the 
PNJL model. We found that D does not have a clear 
order parameter like behavior. However, given the tem- 
perature and the chemical potentials it is possible to esti- 
mate the corresponding D in PNJL model and compare 
with experiments. Such a preliminary comparison has 
been done in this work giving encouraging results and 
indicating the possibility of detecting signatures of ex- 
otic phases in heavy-ion collisions. 
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